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ABSTRACT 

We examine how metallicity affects convection and overshoot in the superadiabatic layer of main 
sequence stars. We present results from a grid of 3D radiation hydrodynamic simulations with four 
metallicities {Z = 0.040, 0.020, 0.010, 0.001), and spanning a range in effective temperature (4950 < 
Toff < 6230). We show that changing the metallicity alters properties of the convective gas dynamics, 
and the structure of the superadiabatic layer and atmosphere. Our grid of simulations show that 
the amount of superadiabaticity, which tracks the transition from efhcient to inefhcient convection, is 
sensitive to changes in metallicity. We find that increasing the metallicity forces the location of the 
transition region to lower densities and pressures, and results in larger mean and turbulent velocities 
throughout the superadiabatic region. We also quantify the degree of convective overshoot in the 
atmosphere, and show that it increases with metallicity as well. 
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1. INTRODUCTION 

The treatment of convection in stellar envelopes is 
one of the largest sources of uncertainty in stellar 
modeling. Convection in stellar models is usually 
describ ed by the Bohm-Vit ense mixing length theory 
(MLT; iBohm-Vitensi [l95l . which_ was derived from 
Prandtl's mixing hypothesis (|Prandt]||1925h . MLT rep- 
resents convection with a single characteristic length, 
which is proportional to the local pressure scale height 
I = aHp, where a is a free parameter. The mixing 
length parameter is arbitrary, and is usually held at a 
constant value, obtained from a calibrated solar model. 
This single parameter formulation of MLT is commonly 
used in ID stellar model calculations, but there are 
other treatments that introduce additional free param- 
eters to account for geomet ric properties of convection 
(see e.g.. lArnett et al.l [20To[ ). While MLT produces an 
accurate st ratification in regions of efficient convection 
(|Chan fc Sofiai , 19891. it does not account for the dynam- 
ical effects of convection, such as turbulent pressure and 
asymmetry in the velocity field, and the approximation 
breaks down near the surface where convection is ineffi- 
cient. 

By fixing the mixing length parameter, the MLT treat- 
ment of convection eliminates dependence on stellar 
properties that are expected to be relevant to convective 
dynamics. These properties include the surface gravity 
and effective temperature of the star, and the chemical 
composition of the co nvection zone. The recent work of 
iBonaca et al.l ()2012[ ) examined how the mixing length 
parameter would need to change to satisfy stellar mass 
and radius constraints from Kepler data. Their results 
indicate that the mixing length parameter changes with 
metallicity. While useful for its application to stellar 
modeling, the results are limited by the unrealistic MLT 
representation of convection. For the past few decades, 
radiation hydrodynamic simulations (RHD) have been 
used to provide a more accurate representation of near 
surface convection. 

Pioneered bv iNordlundl (|1985[ ). 3D RHD simulations 
account for the realistic transition from convective to 
radiative energy transport. Since then, simulations 



have been applied to dwarf stars (e.g. IRamfrez et afl 
[200l IWende et all [200l llYevtag fc Steffen ' '20041^. gi- 



ants (e.g. iCoUet et al.ri2007t iChiavassa et al. 20l3; 
iLudwig fc KucinskasI \20m ~ and several more tar- 
geted studies of individu a l star s fe.g. IRobinson et alj 



200l IStraka et all [2051 \2MA ILudwig et al.l 120091: 



Behara et al.ll2010l ) 

Efforts to systematically study the vari ation of stellar 
convecti on h av e bee n carried out b y 'Ludwig et al.l 
(il99~5L IT9981 [T999I ). iFrevtag et al.l (1999) .^ and 
iTrampedach fc SteinI (|2011D . and there are several 
ambitious efforts to examine metallicity and convection 
using grids of simulations by including a metallicity 
dimension in the grid parameter space. These include 
the ongoing CIFIST Q Ludwig et al. 2009)_and STAG- 
GERGRID (Collet, Magic fc Asplundl 120111) projects. 
Our work follows in a similar vein and employs a grid 
of 3D RHD simulations to isolate the effect of varying 
metallicity while fixing other dimensions of parameter 
space. 

Determining how realistic stellar surface convection is 
sensitive to metallicity remains an area of active research. 
Some of the simulations mentioned above were carried 
out at different metallicities and compositions, and there 
are a few studies com paring simula t ions a t varied metal- 
licity. For example, lJung et al.l (|2007[ ) examined 3D 
RHD simulations of convection in the sun and a pop- 
ulation II star. They find differences in the turbulent 
pressure and kinetic energy, but the stellar parameters 
corresponding to the simulations have different metal- 
licity, effe ctive temperatur e, and surface gravity. Simi- 
larly, Stef fen et al.l ()2009[ ) compared 2D simulations of 
magneto-convection for the Sun and a metal poor solar- 
like analog, but with the same surface gravity and effec- 
tive temperatures. Their preliminary results show that 
the primary effect of the lower opacity in the metal-poor 
star is increased pressure and magnetic field strength. 
We aim to perform a more systematic and rigorous study 
by using a grid of 3D simulations to isolate the effect of 
metallicity over a range in effective temperature. 

Other studies have focused on elemen tal abundances 
and spectral features. For example, lAsplund et al.l 
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(|1999[ ) used convection simulations corresponding to two 
halo stars in conjuction with spectral synthesis to es- 
timate chemical abundances, which indicated that ID 
models have led t o overestima t ed pr imordial Li abun- 
dances. Similarly iCoUet et alJ ([2007) found differences 
in the predicted line strengths from 3D simulations of 
metal-poor red giant stars compared to corresponding 
ID model atmospheres. Their simulations cover a range 
in effective temperature and metallicity, and they found 
that spectral lines of neutral metals and molecules ap- 
pear stronger in 3D simulations than in ID models. 
Results of such studies consistently show cooler atmo- 
sp heres in lowe r meta llicity simulations. A recent study 
by iCollet et al.l ()201lD examined 3D simulations of low- 
metallicity red giants, and confirmed that the cooler at- 
mospheres in low-Z simulations was not the result of the 
approximations to the treatment of scattering. 

Most previous studies of the effect of metallicity on 
convection have focused on atmospheric properties. We 
aim to examine how metallicity determines the structure 
of the convection zone, particularly in the superadiabatic 
layer where energy transport transitions from convective 
to radiative, with the intent of improving stellar mod- 
els. In the present study, we endeavor to systematically 
explore the effect of varying heavy metal abundance and 
convection in the superadiabatic layer by computing sev- 
eral sets of simulations at various metallicities, spanning 
a range in effective temperature. The effective tempera- 
ture range for each metallicity set overlaps, enabling the 
comparison of simulations at varied Z but fixed log(g) 
and log(Tcff). Section [5] outlines the numerical method 
and details of the simulation code, while Section [3] de- 
scribes the individual simulations that we present. In 
Sections |4] through [6] we present the effect of metallicity 
on the thermal structure and gas dynamics. 

2. THE RADIATION HYDRODYNAMICS CODE 

Our simulation code solves the compressible Navier- 
Stokes equations with radiative transfer. Details are pro- 
vided in Tanner et al. (201 ^. It is an updated vers ion of 
the code described in Robinson et al.l ()20 03'. 2004 1), and 
is based on the code of Ki m, Chan & Sofia (jChan h Sofial 
[19891: IKim fc ChanllT998h . 

The simulation domain is a Cartesian box with pe- 
riodic boundary conditions at vertical walls and closed 
surfaces at the top and bottom. Our experiments with 
different boundary conditions consistently show that the 
effect of the boundary layer is less than one scale height 
from the edge of the computational domain, and this 
boundary layer has been removed from the figures. The 
computational domain spans the superadiabatic layer, 
which is the region of transition from convection to ra- 
diation. The bottom of the domain is fully convective, 
while the top is fully radiative. Thus, the computational 
domain spans the transition region from convective to ra- 
diative energy transport, where convection is inefficient 
and MLT is insufficient. 

The extent of the simulation domain is approximately 
9 pressure scale heights (typically less than 0.1% of the 
stellar radius for the simulations presented here), which 
is very small relative to the full extent of the convection 
zone (which is approximately 20 pressure scale heights in 
the sun) , so curvature and radial variation of gravity are 
negligible and can be ignored. Although the simulation 



domain is small relative to the full extent of the con- 
vection zone, it encompasses precisely the region where 
MLT fails to capture the effect realistic convection, such 
as the additional support of turbulent pressure. 

Radiative transfer is treated with the diffusion approx- 
imation in the optically thick part of the domain. In the 
optically thin layers where the diffusion approximation 
is invalid, the code can solve for the m ean intensity us- 
ing t he 3D Eddington approximation (jUnno fc Spiegell 
Il966f ) or integration over long ray characteristics (simi- 
lar to the method described by Stein & Nordlund ( 2003) . 
The simulations were computed using the 3D Eddington 
solver, which compares well with long characteristic ray 
integration, and offers some computational performance 
benefits. A detailed comparis on of the radiative tr ansfer 
solvers has been presented in lTanner et al.l (|2012D . The 
code uses the OPAL equation of state and opacity ta- 
bles ( Rogers fc Navfonovl 120021 ) and the IFerguson et al.l 
(|2005l ) opacitv tables at low temperatures. 

3. GRID OF SIMULATIONS 

Each simulation is characterized by its energy flux, sur- 
face gravity, and composition. It is important to note 
that while we directly set the surface gravity and com- 
position, the energy flux is a quantity that is calculated. 
We define the effective temperature of the simulation as 

^cff - — (1) 

where cr is Stcfan-Boltzmann constant, and the radiative 
flux, i^iad, is determined by the radiative transfer solver. 
Changing the metallicity in a simulation alters the ther- 
mal balance in the atmosphere, resulting in a different 
energy flux. Without precise control over the radiative 
energy flux, it is difficult to compare simulations at con- 
stant effective temperature and varied metallicity. To 
overcome this, we have computed a grid of simulations 
for several metallicities that span a small range in effec- 
tive temperature. 

The grid comprises four sets of simulations, each at a 
different metallicity. Each set spans a small range in ef- 
fective temperature, which overlap with each other. All 
simulations in the grid have the same surface gravity set 
to log (g) = 4.30. Properties of the grid are listed in Ta- 
ble [TJ We can isolate the effect of metallicity by interpo- 
lating within the grid to achieve a desired effective tem- 
perature, or by comparing simulations that happen to 
have similar energy fluxes. Simulations use opacity and 
equation of state tables with the heavy e lement abun- 
dance mixture of iGrevesse fc Sauvall (|1998[ ). 

Since the computational domain spans a sufficiently 
small fraction of the stellar radius so that the gravita- 
tional field is homogeneous and curvature is neglected, 
the simulations do not contain information about the 
stellar radius. As a result, the zero point of the radial 
coordinate is arbitrary, and the simulation domain can- 
not be expressed as a fractional stellar radius. Instead 
of using radius, we typically present properties of the 
simulations relative to the effective temperature surface. 
This reference point (where (T) = Toff) can be inter- 
preted as the stellar 'surface' because it is the point at 
which a black body yields the same energy flux as that of 
the simulation. We define the height such that the zero 
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TABLE 1 

Properties of the simulations in the grid. All simulations 

HAVE THE SAME SURFACE GRAVITY (logg = 4.30). 



T D 
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Zi 
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-'eft 


log eft 


sl 


0.040 


0.715 


4947 


3.694 


s2 


0.040 


0.715 


5201 


3.716 


s3 


0.040 


0.715 


5461 


3.737 


s4 


0.040 


0.715 


5719 


3.757 


s5 


0.020 


0.735 


5121 


3.709 


s6 


0.020 


0.735 


5370 


3.730 


s7 


0.020 


0.735 


5626 


3.750 


s8 


0.020 


0.735 


5897 


3.770 


s9 


0.010 


0.745 


5326 


3.726 


slO 


0.010 


0.745 


5569 


3.746 


sll 


0.010 


0.745 


5803 


3.764 


sl2 


0.010 


0.745 


6027 


3.780 


sis 


0.001 


0.754 


5735 


3.759 


sl4 


0.001 


0.754 


5896 


3.771 


sl5 


0.001 


0.754 


6063 


3.783 


sl6 


0.001 


0.754 


6235 


3.795 



point is fixed at the Toff surface, and it increases toward 
smaller optical depth. 

In the following sections we compare space- and time- 
averaged quantities over the vertical (radial) coordinate. 
The horizontal and temporal mean over a time At and a 



box cross-sectional area of L^r. x Ly is: 



(9) 



1 1 

Ai LxLy 



qdxdydt. (2) 



Note that quantities are averaged over a sufficient in- 
terval to obtain statistical convergence, and the start of 
the time interval is after the simulation has thermally 
relaxed. The rms of the same quantity is: 



(3) 



The correlation function between two turbulent quan- 
tities qi and q2 is defined as: 



C[qi,q2] = 



(gig2) - (gi)(g2) 

'?l,rms92,rms 



(4) 



4. RESULTS: SAL STRUCTURE 
4.1. Mean Stratification 

The structure of the near surface layers is determined 
largely by the inefficient convection. It is in precisely this 
region that convection models such as MLT fail to pro- 
vide accurate stratifications. Figure [T] shows the density 
as a function of total pressure (which includes both gas 
and turbulent pressure) for four simulations. 

The simulations in Figure [T] were selected to show the 
effect of increasing the energy flux and changing the 
metallicity. The region of transition from convective to 
radiative energy transport, known as the superadiabatic 
layer (SAL), occurs just below the stellar surface. In 
Figure [1] the SAL is apparent as an abrupt change in the 
density-pressure relationship. Over this pressure range, 
increasing the effective temperature primarily changes 
the density in these layers, while changing metallicity 
substantially changes both the density and the location 



of the SAL. Higher metallicity simulations have SALs 
that occur at lower pressure and density. 

The density in the SAL reflects the change in opacity 
as a result of changing metallicity. At these tempera- 
tures, the primary opacity source is from the H~ ion. In- 
creasing the metallicity results in low ionization elements 
donating more electrons to hydrogen atoms, thereby in- 
creasing the opacity. The energy flux is determined by 
pn and the source function, so to maintain the same en- 
ergy flux at higher opacity (k), the density (p) must be 
correspondingly lower. 
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Fig. 1. — Structures through the SAL for four select simulations. 
Simulation pair (sl3,sl4) shows how the stratification changes 
with effective temperature, while simulations pairs (s4,sl3) and 
(s8,sl4) have essentially the same effective temperatures but differ- 
ent metallicities. The density and pressure of the SAL are sensitive 
to changes in metallicity. 



4.2. Superadiabatic Excess 

Near the surface, where convection is inefficient, a use- 
ful quantity to examine is the radial variation of the su- 
peradiabatic excess (V — Vad), which is the difference 
between the temperature gradient, defined as: 



dln(r) 
dln(P)' 



(5) 



and the adiabatic temperature gradient. Turbulent gas 
dynamics in the SAL contribute a significant fraction to 
the pressure term in Equation [5j The total pressure, 
P = Pgas -l- fturb, IS the sum of the gas and turbulent 
pressures, where the turbulent pressure is calculated with 
the defintion in Section [521 In the deep layers where the 
structure is nearly isentropic and convection is efficient, 
the temperature gradient is essentially adiabatic. Nearer 
to the stellar surface, where the energy flux transitions 
from convective to radiative, we see a departure from adi- 
abaticity characterized by a sharp peak in the superadia- 
batic gradient. The difference between the two gradients 
(i.e. the superadiabatic excess) , provides a measure of the 
inefficiency of convection, and reveals the depth at which 
there is a significant departure from efficient convection. 

Figure [2] shows the superadiabatic gradients as a func- 
tion of pressure for our grid of simulations. Each panel in 
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the figure corresponds to a different metallicity set. The 
temperature structures show a mild systematic change in 
the maximum suberadiabaticity (maximum of V — Vad) 
with the effective temperature of the simulations, al- 
though the variation is markedly less apparent than what 
is predicted by stellar models computed with MLT. The 
more constrained range in superadiabaticity relative to 
M LT is generally con s istent with the grid of simulations 
of iTrampedach et alj ()2010D . which show a smaller vari- 
ation in the superadiabatic peak compared with MLT 
models at evolutionary stages from the main sequence 
to giants. The simulations with larger effective temper- 
atures have higher SALs, particularly for high-Z simula- 
tions, but the maximum suberadiabaticity becomes less 
sensitive to Teg at near-Solar and lower-Z. We note, how- 
ever, that while the temperature ranges overlap, they are 
not identical for each metallicity simulation set. 

When the metallicity is changed and Teff is fixed, we 
find no significant change in the maximum superadia- 
baticity. The top panel of Figure [3] compares the supera- 
diabatic gradients in two simulations that have almost 
the same effective temperatures (the effective tempera- 
tures differ by only 27K) but a drastic change in metal- 
licity (Z differs by a factor of 40). In both simulations 
the maximum superadiabatic excess is approximately the 
same. We do, however, see a marked shift in the depth of 
the SAL. At lower metallicity, the conversion from effi- 
cient to inefficient convection happens at higher pressures 
and densities. 

The degree of superadiabaticity corresponds to the rate 
at which energy transport transitions from convective to 
radiative. A large value of the maximum SAL is the re- 
sult of a steep entropy gradient. The relatively constant 
SAL maximum as a function of metallicity suggests that 
the thermal structure changes (examined in section |4TT|) . 
rather than a large change in the convection-to-radiation 
transition. 

We summarize the simulated SALs in the lower panels 
of Figure |3l which show the maximum superadiabaticity 
and its location for all the simulations in the grid. In 
the top panel, each metallicity set shows a shallow trend 
with effective temperature, but there is no remarkable 
offset between them. Although the maximum value of 
the SAL remains unchanged, the bottom panel of Figure 
[3] shows that the change in metallicity introduces a shift 
in the location of the SAL. Lower-Z simulations have 
correspondingly higher densities to maintain the same 
energy flux. Thus, radiative energy transport starts to 
become significant at higher density in the lower-Z sim- 
ulations. 

5. RESULTS: DYNAMICS OF CONVECTION 

5.1. Mean Convective Velocities 

One of the distinctions between mixing length theory 
and realistic convection is the asymmetry of the veloc- 
ity field resulting from upflows being characterized by a 
larger fllling factor at almost all depths in the simulation 
domain. The hot upflowing gas rises nearl y adiabatically 
at approximately constant entropy (e.g. iLudwig et al.l 
1199a IStein et al.l 119991) until it radiates away energy 
near the surface before cycling back into the convection 
zone. Stellar convection is characterized by large up- 
flowing regions (granules) separated by cool fast down- 



drafts (intergranular lanes). Figure 2] shows the vertical 
velocity at a horizontal slice where (T) — T^s for two 
simulations with different metallicities. The flow pat- 
tern qualitatively resembles solar granulation with large 
expanding updrafts separated by cool fast downdrafts. 
The difference in metallicity, however, causes the granu- 
lation pattern to be visibly different. Granulation in the 
higher-Z simulation shows larger granules and sharper 
inter-granular lanes. 

Asymmetry in the velocity field means that the aver- 
age vertical velocity, (w), is dominated by the upflows 
through most of the simulation domain. As shown in 
Figure [51 the average vertical velocities are small deep in 
the convection zone and rise to a maximum just below 
the surface near the peak of the SAL. The velocity proflle 
is inverted above the convectively stable region (where 
V < Vad), and the average velocity becomes negative. 

The maximum average vertical velocities of the differ- 
ent simulations are compared in FigurelHl The trend with 
effective temperature (from velocity proflles in Figure [5]) 
is apparent, but there is also a distinct metallicity depen- 
dence. At a given effective temperature, simulations with 
higher metallicity exhibit larger convective velocities. 

The trend is a direct result of the changes in density 
stratiflcation, which change in response to the different 
opacities caused by changing the metallicity (see Section 
14.11) . In simulations with the same Tes, convection is 
responsible for transferring the same energy flux. The 
different stratifications, however, require different veloc- 
ities. High-Z simulations have higher opacity, and con- 
sequently lower density through the SAL (see Figure [T]) , 
so the convective velocities must be larger to sustain the 
same convective flux. 

Larger convective velocities in the SAL mean that up- 
flows will remain coherent over a larger spatial scale, even 
in the radiative domain where convection is no longer 
driven. In Section r5.3l we examine the strength of convec- 
tion beyond the convectively stable boundary, and pro- 
vide a measure of how it changes with metallicity. 

5.2. Turbulent Pressure 

One of the physical phenomena absent in MLT-like pre- 
scriptions for convection is the additional pressure sup- 
port from turbulence. The so-called turbulent pressure 
is caused by turbulent fluctuations in the velocity flcld, 
and is a non-negligible contribution to hydrostatic bal- 
ance through the SAL, and causes structural differences 
where turbulent pressure is signiflcant. 

We estimate turbulent pressure from the rms (Equa- 
tion [3]) of the vertical velocity field as: 

-Pturb = PWrms- (6) 

Simulations of solar convection estimate turbulent 
pressure support at approximately 15% t o 18% of the 
gas pressure near the peak of the SAL (|Kupkal I200I 
IBeeck et al.l IMl) . We see similar results from our 
code. Turbulent pressure is negligible deep in the near- 
adiabatic part of the convection zone because of ex- 
tremely large gas pressure, and is also quite small in the 
optically thin region where the gas density and rms ve- 
locities are small. 

Figure [7] compares the radial variation of turbulent 
pressure (relative to gas pressure) from all of the sim- 
ulations. Turbulent pressure is a few percent in the deep 
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Fig. 2. — Superadiabatic gradients (V — Vad) for the grid of simulations. Each panel shows a set of simulations over a small range in 
effective temperature and at a given metaUicity. At a given metaUicity, the maximum superadiabaticity changes with T^ff , but at lower 
metallicity the maximum superadiabaticity becomes less sensitive to T^g. 

mixing. The amount of additional mixing, or overshoot, 
win depend on the stratification of the atmosphere and 
the convective dynamics. 

Simulations treat convection and the effect of over- 
shoot self-consistently, and make no distinction between 
them. We can, however, quantify the amount of over- 
shoot by measuring correlations between turbulent fluc- 
tuations. In convectively unstable regions, the upflows 
are coherent and turbulent fluctuations are highly corre- 
lated. By measuring the vertical (radial) location where 
the fluctuations are no longer correlated, we can deter- 
mine the height at which the convection no longer resem- 
bles the behavior of the unstable region. 

We define the end of the convectively unstable region 
(and the start of the overshoot region) to be the height at 
which V — Vad is zero (see Figure [2]). We define the end 
of the overshoot region to be where fiuctuations in tem- 
perature and vertical velocity are no longer correlated, 
where the degree of correlation is defined by Equation 
H This definit ion ha s been used bv lChan fc So fia (I989|) 
and [Robinson et al.l ()2004[) . The size of the overshoot 
region is the difference between these two heights. 

The correlation through the deeper part of the convec- 
tion zone is quite constant over this range in T^s and Z, 
ranging from C[T, Vz] = 0.77 — 0.85, which is similar to 
the value measured by Chan & Sofia (1989). The rate 
at which the correlation drops, however, is sensitive to 
metallicity. Figure |9] compares the temperature- velocity 
correlations and superadiabatic gradients in two simu- 
lations with different metallicity and the same effective 



convection zone, and increases at smaller depth, reaching 
a maximum near the peak of the SAL. Turbulent pres- 
sures in our grid can be as high as 10-16% of the gas 
pressure. 

As with the mean vertical velocities, we find that the 
turbulent pressure support scales with both the energy 
flux and metallicity. The high-Z simulations show more 
support from turbulence (as a fraction of gas pressure) 
relative to the low-Z simulations. Comparing the tur- 
bulent pressure from each simulation (Figure 18]) shows 
that the maximum turbulent pressure visibly increases 
with effective temperature, and there is a clear system- 
atic offset that is correlated with changes in metallicity. 
The increase in turbulent pressure is a result of larger 
rms velocities, which are also correlated with effective 
temperature and metallicity. 

5.3. Overshoot 

Convective regions in ID stellar models are well defined 
and unambiguous. The radial domain that is unstable to 
convection can be defined according to many equivalent 
criteria, such as Vad < V, ds/dz < 0, or Fconv > 0. 
Stratifications from 3D RHD simulations do not provide 
such clear boundaries, and there are several ways to esti- 
mate the edge of a driven convection zone, or the extent 
of a region influenced by convection. 

It is well known that updrafts would not simply stop at 
the edge of the convectively unstable layer. Instead, their 
momentum would carry them some distance beyond the 
convective boundary effectively extending the region of 
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Fig. 3. — (a): Superadiabatic gradients from simulations s4 and 
sl3, which have different metaUicities but the same effective tem- 
perature. The maximum of the SAL remains essentiaUy unchanged 
despite the considerable difference in metallicity. The location of 
the SAL is sensitive to metallicity. (6): The maximum value of 
the superadiabaticity and (c) location of the maximum supcradi- 
abaticity for all simulations in the grid. Each line in panels (b) 
and (c) corresponds to a set of simulations with a fixed metallicity. 
Simulations s4 and sl3 are marked with diamond symbols. Metal- 
licity does not appear to directly affect the maximum of the SAL, 
but it does introduce a shift in the location of the SAL with respect 
to density. 

temperature. The zero point of the radial height is de- 
fined to match the average effective temperature surface. 
The edge of the convectively unstable region, as defined 
by where the superadiabatic gradient is zero, is almost 
unchanged in the two simulations. The region affected by 
overshooting, as measured by the temperature-velocity 
correlation, is clearly larger in the higher metallicity sim- 
ulation. 

Figure [TOl shows the size of the overshoot region for the 
entire set of simulations. There appears to be a small 
correlation of overshoot with effective temperature, par- 
ticularly in the high-Z simulations, but the dominant ef- 
fect is clearly by varying metallicity. The larger degree 



of overshoot in high-Z simulations is attributed to the 
larger mean and turbulent velocities. The faster moving 
high-Z upflows carry more momentum and are able to 
traverse more space before losing their identity. 

6. RESULTS: T-T RELATIONS 

Convection also affects the structure above the region 
of driven convection because of convective overshoot. By 
directly solving the hydrodynamics equations, simula- 
tions provide a self-consistent representation of the effect 
of convection, including overshoot, on the atmosphere. 
The primary way in which the stratification in and above 
the convective envelope is altered is from the additional 
sup port provided by turbulent pressure (described in Sec- 
tion [5?2]). 

ID Stellar structure models set the outer bound- 
ary condition by imposing a structure defined by a 
T-T relation. The Eddington T-r is often used, but 
semi-empirical relations, such as those provi ded by 
Krishna Swa^ (|1966D and lVernazza et al.l (|198lD (VAL) 
can also be applied. The semi-empirical relations implic- 
itly include effects such as turbulent pressure, but are 
only valid for the Sun. The location of the T^s surface 
with respect to optical depth is defined by the T-r rela- 
tion. The Eddington T-r defines the Tea surface to be 
at r = 2/3, while the semi-empirical relations show the 
Toff surface at smaller optical depths. Examining T-r 
relations from simulations demonstrates how convective 
turbulence and overshoot affects the atmospheric struc- 
ture, and to what degree it is sensitive to changes in 
metallicity. 

The variation in T-r relations with energy flux is rather 
small over the range of effective temperatures here. In 
simulation sets of a given metallicity, the temperature 
variation at a given optical depth is typically less than 
a few percent. The T-r relation is more sensitive to 
changes in metallicity, however, as demonstrated in Fig- 
ure 111! which shows T-r relations for two simulations 
with the effective temperature but different metaUicities. 
The Eddington, Krisnha Swamy, and VAL T-r relations 
are included for reference. It is immediately apparent 
that the effect convective overshoot is sensitive to metal- 
licity. The comparison shows a discrepancy in the tem- 
perature at a given optical depth of up to 6%. 

There is a clear systematic effect in response to the 
changing metalicities. Low-Z stratifications are cooler 
than high-Z stratifications at a given optical depth. The 
example presented in Figure [TT] shows that the high-Z 
and low-Z atmospheres have temperatures of approxi- 
mately T/Tcff = 0.846 and 0.793 at r = 0.01, respec- 
tively. 

Cooler atmospheric tempe ratures at low-Z in 3D sim- 
ulations have been seen by lAsplund et all (11999) and 
iCollet et al.l (j2011| ) . They explain that the thermal struc- 
ture of the atmosphere is determined by balancing the 
competing effects of radiative heating and adiabatic cool- 
ing from expanding upflows. The decreased metallicity 
simulation has lower opacity and is less able to absorb 
radiation. This leads to less radiative heating, so that 
balance is achieved at lower temperatures. 

7. DISCUSSION 

Simulations of near-surface convection provide a realis- 
tic description of convection, and self-consistently couple 
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Fig. 4. — Vertical velocity at the effective temperature surface showing the characteristic granulation pattern from snapshots 
from two simulations in the grid. The two simulations have the same surface gravity and effective temperature, but metallicities 
of Z = 0.040 (s4, left) and Z = 0.001 (sl3, right). 
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Fig. 5. — Radial variation of average vertical velocities for all simulations in the grid. Velocities are small in the deep 
region, and reach a maximum near the peak of the superadiabaticity. Larger velocities are measured in simulations with 
higher metallicity or effective temperature. 
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Fig. 6. — {top): Radial variation of vertical velocity for simula- 
tions s4 and sl3 (identified with diamond symbols in the bottom 
panel), which have approximately the same T^ff. {bottom): The 
maximum average vertical velocity for all simulations in the grid. 
The maximum velocity increases with effective temperature and 
metallicity. 

the atmosphere and superadiabatic layer to the adiabatic 
region below. The resulting stratification differs from, 
and is superior to ID stellar models computed with MLT 
treatment of convection. Differences between simulations 
and ID models arise because the MLT treatment of con- 
vection lacks a description for some physical processes 
that are present in the simulations. In particular, MLT 
does not include the turbulent contribution to pressure 
support. Additionally, the upflowing granulation near 
the surface modifies the T-t relation as a result of en- 
hanced adiabatic cooling that ID models do not include. 
Providing a representation of processes like these in ID 
stellar models could potentially improve their accuracy. 

In Sections H] and [S] we have shown how metallicity af- 
fects near-surface convection in stars. We are able to 
isolate the effect of metallicity by using a grid of simu- 
lations comprising sets of simulations at a fixed Z, but 
with overlapping Teg ranges. Increasing metallicity in- 
creases the opacity, which alters the way in which energy 
transport transitions from convective to radiative. 

When the effective temperature is held fixed, changing 
the metallicity alters the mean stratification through the 
superadiabaitc region, with low-Z simulations exhibiting 
higher densities and pressures through the SAL. Addi- 
tional changes in pressure come from the contribution 
to hydrostatic equilibrium from turbulent pressure. Be- 
cause of the larger gas pressure at low-Z, the turbulent 
pressure relative to the total pressure becomes smaller. 



although it is still significant. The increase in density is 
required to compensate for the drop in opacity, in order 
to maintain the radiative flux. This effect is expected, 
and certainly not specific to 3D simulations, but it has 
ramifications for other properties of convection present 
in the simulations. 

The simulations show an adjustment of the superadi- 
abatic excess as a result of changes to the mean strat- 
ification. In the deep part of the convection zone, en- 
ergy is transported almost entirely by convection and the 
temperature gradient is essentially adiabatic. Nearer to 
the surface, the temperature gradient becomes superadi- 
abatic as energy is transferred outward more efficiently 
by radiation. The degree of superadiabatic excess marks 
the transition from convective to radiative energy trans- 
port, and is sensitive to changes in metallicity. 

Varying metallicity in the simulations induces remark- 
ably little change in the maximum value of the supera- 
diabaticity, even when the metallicity is reduced by a 
factor of 40. The extent of the superadiabatic region 
also shows little remarkable variation between simula- 
tions. This suggests that the rate of transition from con- 
vective to radiative energy transport is not drastically 
different in simulations of different metallicities. The lo- 
cation of the transition, however, is quite different, as 
evidenced by a shift in the location of the superadiabatic 
peak. The SALs in the low-Z simulations occur at higher 
density. The low-Z simulations have smaller opacity at a 
given density, and as a result the start of the transition 
to radiative energy transport occurs at higher density. 

The structure at small optical depth is set by the bal- 
ance of radiative heating and adiabatic cooling from the 
convective updrafts. 3D simulations are able to account 
for accurate adiabatic cooling from the asymmetric fiow 
that his characteristic of granulation. Lower opacity in 
low-Z simulations results in less radiative cooling, and 
causes steeper T-t gradients and cooler atmospheres. 
This effect has been previously demo nstrat ed using 3D 
simulations by lAsplund et al.l ()1999D and iCollet et al.l 
(2007, 20l3). When metallicity is changed, our simu- 
lations show a substantial deviation (Figure [TT] shows a 
change in temperature of up to 8%) from the simple T-t 
relations commonly used in stellar modeling. 

Convective velocities are determined by the stratifi- 
cation, and the amount of energy flux that is trans- 
ported. Simulations with higher effective temperatures 
must transport more energy flux, and the convective ve- 
locities are expectedly larger. The correlation of velocity 
and Teff exhibited in our simulation is consisten t with 
the trend reported by iTrampedach fc SteinI (|201lD . The 
energy fluxes are the same in simulations with the same 
Toff (but different Z), so differences in the convective ve- 
locities is a result of differences in the stratiflcations. Be- 
cause the high-Z simulations have lower density, the con- 
vective velocities must increase to maintain the energy 
flux. Turbulent velocities are also strongly correlated 
with metallicity. Our simulations show that the turbu- 
lent pressure as a fraction of gas pressure, which is de- 
termined by the density and turbulent (rms) velocities, 
increases with effective temperature and metallicity. 

Most of the simulation domain is convectively unstable, 
and, furthermore upflowing plumes do not simply stop 
in the radiative layers. Instead, the momentum of the 
upflows carries them a certain 'overshoot' distance be- 



9 




Height (km) Height (km) 

Fig. 7. — Radial variation of average turbulent pressure for all simulations in the grid. Turbulent pressure follows the same 
trend as the mean velocities (Figurej^l, increasing with effective temperature and metallicity. 



fore cycling back into the convection zone. An expected 
consequence of larger mean and rms velocities is an in- 
crease in the overshoot distance. We are able to measure 
boundaries on the overshoot region using correlations be- 
tween turbulent fluctuations in temperature and velocity. 
We find that high-Z simulations show significantly more 
overshoot. Our simulations show a consistent and clear 
correlation between metallicity and overshoot, and in the 
most extreme case (between simulations with Z = 0.040 
and Z — 0.010) we find a factor of 4.5 change in over- 
shoot distance. 

The effect of convection, and its variation with stellar 
evolution, is typically not accurately represented in stel- 
lar modeling. Most stellar models are calculated using 
MLT, or a similar prescription for convection, but there 
have been several attempts at includi ng realistic con- 
vectio n in stellar models. For example, iDemargue et all 
(|1997f) introduced a variable mixi ng length parame - 
ter derive d from the simulations of iKim et al.l (|1996( ). 
[Rosenthal et al\ ([1999) matched an envelope model to a 
solar simulation to include the effect o f turbulent con - 
vection on oscillation frequencies, while iLi et al\ ()2002f ) 
included a two-parameter model for turbulent pressure 
and kinetic energy (calibrated with simulations) to ID 



stellar models. 

The goal of simulating stellar envelope convection re- 
mains a generalized improvement over the current MLT- 
like treatments. Our simulations show that that im- 
posing a solar-calibrated mixing length parameter on all 
stellar models is not accurate, and ultimately, any gen- 
eralized improvement for the treatment of convection in 
stellar models will have to account for composition. Sim- 
ulations can help improve stellar models by providing 
a self-consistent and realistic description of convection 
across the relative parameter space, which includes com- 
position. While metallicity is an important contribution 
to the composition, other significant factors, such as he- 
lium abundance, warrant consideration as well. To that 
end, we are currently using simulations to investigate the 
effect of helium on convection. 

JT is supported by NASA ATFP grant 
#NNX09AJ53G to SB, and acknowledges a PGS-D 
scholarship from the Natural Sciences and Engineering 
Research Council of Canada. This work was supported 
in part by the facilities and staff of the Yale Univer- 
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Fig. 9. — Convective overshoot for simulations s4 (top) and sl3 (bottom). The left side shows the superadiabatic excess and 
correlation functions that are used to define the region of overshoot. The right side shows the slices of the vertical velocity 
near the T^g surface. Dashed lines show the measured boundaries of the overshoot region. The high-Z simulation (s4) has 
larger velocities and more overshoot. 




Fig. 10. — The extent of convective overshoot for all simulations 
in the grid. High-Z simulations have larger velocities which result 
in more convective overshoot. 



12 



1 — I — I — I I 1 1 1 — '■ 

VAL Model "C" 
Krishna Swamy 
Eddington 

■ Z=O.OOI log(T,ff) = 3.757 

■ Z=0.040 log(T„ff) = 3.759 



/ 




7 8 9 

0.0 1 



7 8 9 

O.I 



1 



Optical Depth 

Fig. 11. — Average T-t relations for simulations s4 and sl3 with 
varied metallicity, compared with analytic and semi-empirical T-t 
relations for the Sun. Low metallicity simulations have less radia- 
tive heating, resulting in a steeper T-t relation in the optically 
thin layers. Dotted lines mark the T^ff surface and the r = 2/3 
optical surface, which by definition, occur at the same depth in the 
Eddington T-t relation. 



